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NOMENCLATURE 


a 

A 


D 

F P 

FS 

F V 

k 

L 

m 

M 

Nu 

P 

q 

Q 

R 

Re 

RTR 

S 


speed of sound, m/sec 
cross-sectional area of sensor wire, m 2 
specific heat, J/kg °K 


circumference of sensor wire, m 


pressure coefficient, 


P L ~ P oo 


, dimensionless 


diameter of sensor wire, m 

view factor for planetary radiation, dimensionless 
view factor for solar radiation, dimensionless 
view factor for vehicle-emitted radiation, dimensionless 
convective heat-transfer coefficient, W/m 2 °K 
thermal conductivity, W/m °K 
length of sensor wire, m 

parameter defined by equation (15), per m 

V 

Mach number, — , dimensionless 

a 

Nusselt number, , dimensionless 
h-a 

pressure, N/m 2 

pV 2 

dynamic pressure, — N/m 2 

heat-transfer rate, W/m 2 

gas constant for a particular gas, J/kg °K 


Reynolds number, 


pVD 


, dimensionless 


77? 

recovery temperature ratio, — , dimensionless 

T t 

solar constant, W/m 2 



ratio of length measured along forebody surface, from axis of vehicle, to maximum 
radius of forebody 

time, sec 

temperature, °K 

speed, m/sec 

distance along sensor wire, m 

absorptance of sensor wire for radiation emitted by a particular source, dimensionless 

c 

ratio of specific heats, dimensionless 

c v 

emittance of a particular source of radiant energy, dimensionless 
dynamic viscosity, N sec/m 2 
density, kg/m 3 

Stefan-Boltzmann constant, W/m 2 °K 4 
time constant, sec 

Subscripts and Special Notation 

atmosphere 

convection 

local 

planet 

recovery 

sensor 

Sun 

total 


vehicle 



(x,t) 

OO 


location x on sensor wire at time t 
free stream 


ANALYTICAL DESIGN OF SENSORS FOR MEASURING ATMOSPHERIC TEMPERATURE 


DURING TERMINAL PHASE OF PLANETARY ENTRY 

John P. Millard, Michael J. Green, and Simon C. Sommer 
Ames Research Center 


SUMMARY 


An analytical study was conducted to develop a sensor for measuring the temperature of a 
planetary atmosphere from an entry vehicle traveling at supersonic speeds and having a detached 
shock. Such a sensor has been used in the Planetary Atmosphere Experiments Test Probe (PAET) 
mission and is planned for the Viking-Mars mission. The study specifically considered butt-welded 
thermocouple sensors stretched between two support posts; however, the factors considered are 
sufficiently general to apply to other sensors as well. This study included: (1) an investigation of the 
relation between sensor-measured temperature and free-stream conditions; (2) an evaluation of the 
effects of the extraneous sources of heat; (3) the development of a computer program for evaluating 
sensor response during entry; and (4) a parametric study of sensor design characteristics. 

The results of the parametric study show that inaccurate knowledge of the recovery 
temperature ratio is the potential source of the greatest error for the sensor designed for the 
PAET-Earth mission with deployment at Mach number 2. The effects of radiation are negligible for 
this particular mission. For deployment at Mach number 3, the sensor does not track recovery 
temperature as well because the temperatures are high and changing rapidly with time. If the sensor 
wire designed for PAET were shortened by a factor of 2, errors induced by conduction effects 
would increase by about 2° K. However, if this short sensor were used for the Viking-Mars mission, 
errors of the order of 50° K would result. Thus, it cannot be assumed that a sensor designed for one 
mission will be adequate for another. 


INTRODUCTION 


An atmospheric temperature sensor was selected as one of the instruments to be flown aboard 
the Planetary Atmosphere Experiments Test (PAET) vehicle. This vehicle was to enter the Earth’s 
atmosphere at a velocity of 7000 m/sec with the mission of evaluating techniques for measuring 
atmospheric composition and structure. The temperature sensor was to provide a measurement of 
free-stream ambient temperature, as a function of altitude in the terminal phases of entry, to an 
accuracy of 1 0° K or better. The sensor was to be deployed when the vehicle had slowed, because 
of drag to about Mach number 2 or 3, after which it was to remain operational until impact. The 
sensor was positioned on the front face, near the base of the vehicle forebody, and extended 
beyond the boundary layer but behind the main shock wave. 



The sensor selected was a butt-welded thermocouple stretched between two support posts. Its 
advantages are (l)the wire diameter can be specified to yield a tailormade time response for a 
specific mission (wires as small as 0.025-mm dia could be used), (2) the temperature at the 
butt-weld junction can be made independent of support-post temperature by the use of sufficiently 
long wires, (3) metals usually have a higher reflectance for solar and infrared radiation than do other 
materials, thus reducing solar and infrared inputs; and (4) the cylindrical configuration of a rolled 
butt-welded joint is a configuration for which there is much heat-transfer data in the literature. 

An analytical study was undertaken as a basis for the design of a butt-welded thermocouple 
sensor that would yield reliable data. This study included: (l)an investigation of the relation 
between sensor-measured temperature and free-stream conditions; (2) an evaluation of the effects of 
extraneous sources of heat; (3) development of a computer program for evaluating sensor response 
during entry; and (4) a parametric study of sensor design characteristics. 

This paper documents the analytical study and illustrates its value. The factors considered and 
the computer program used are sufficiently general to apply to similar missions. Numerical 
examples pertain to possible Earth or Mars missions; PAET was to be tested first in the Earth’s 
atmosphere, and a similar temperature sensor is planned for the Viking-Mars mission. 


FACTORS INFLUENCING MEASUREMENT 


The measurement of free-stream ambient temperature with a sensor located behind a shock 
wave is a complex problem. The ambient temperature must be determined indirectly from a 
measurement of total temperature, which for reference is defined in equation (1) for a gas with 
constant ratio of specific heat 7. In the absence of radiation, total temperature is constant across a 
shock, and therefore the value at the sensor location is the same as in the free stream. Density and 
velocity vary with position behind a shock, however, and these influence sensor response time and 

efficiency of conversion of kinetic energy 
of the gas to thermal energy at the sensor. 
Thus, it is necessary to know the flow 
characteristics at the sensor location. 

Another problem is that of extraneous 
heat sources, which can render a measure- 
ment meaningless: the Sun, planetary and 
vehicle infrared radiation, infrared emission 
of the sensor itself, and conduction 
through support leads. The importance of 
these heat sources can be evaluated only by 
comparing their heat inputs with that from 
the air stream, and their implications for 
sensor performance will vary with the 
particular mission. 



Figure 1. — Temperature sensor designed for PAET. 


In this section, factors that influence 
measurement are discussed in some detail 
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and incorporated into an energy equation 
for relating free-stream ambient tempera- 
ture to sensor response. This discussion is 
general, but numerical values are presented 
that pertain to butt-welded thermocouple 
sensors and Earth and Mars missions. The 
sensor designed and constructed for PAET 
is shown in figure 1 ; the direction of flow is 
perpendicular to the plane of the page. For 
redundancy, this particular sensor contains 
three thermocouples mounted in parallel; 
for analysis, only one wire need be con- 
sidered. Figure 2 illustrates the blunt-body 
PAET vehicle and its detached shock wave. 
This is a shadowgraph of a PAET model in 
a ballistic range at a Mach number of 3. 1 
Tables 1 and 2 give trajectory and plane- 
tary data for the PAET and Mars missions. 


TEMPERATURE 


/ SENSOR 
/ LOCATION 



Figure 2. - Blunt-body PAET vehicle and detached 
shock wave. 


Total Temperature 


Total temperature is the measurement from which free-stream ambient temperature can be 
deduced. The relation is 

T=T 1 1 M 2 ] (1) 

t 00 \ 2 °° J 

The term on the right is the sum of ambient temperature and the temperature increment contributed 


by kinetic energy of the gas molecules 
relative to the vehicle. An important conse- 
quence of this equation is that a vehicle 
undergoing rapid deceleration from high 
Mach numbers will encounter large and 
rapid changes in total temperature. This is 
illustrated in figure 3 for Earth and Mars 
missions. The measurement problem in- 
cludes the accurate tracking of this rapidly 
varying total temperature (41° K/sec at 
M = 2 for PAET). 

In the design of the sensor, the free- 
stream ambient temperature profile and 
vehicle Mach number are assumed to be 
known; the total temperature versus time is 



Figure 3. — Total temperature versus time for PAET 
and Viking missions. 


Provided by courtesy of Mr. Robert I. Sammonds, Hypersonic Free Flight Branch, Ames Research Center, NASA. 
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computed from equation ( 1 ), and the response of a candidate sensor to that imposed temperature 
variation is sought. 


Flow Conditions at Sensor Location 


Values of Mach and Reynolds numbers at the location of the sensor influence the speed and 
efficiency of sensor response. In the following, variables necessary to evaluate these numbers are 
defined in terms of free-stream conditions. The definitions were obtained or deduced from 
reference 1 . 


Local Mach number can be written; 


M L = 


where 


For Moo < 1 


and foryWoo > 1 


i 

Pl.oo 

7 




1/2 


P Jl 

Pt, 


,l_\ ( 7+1 )Moo 2 *1,7-1 r 7+1 -| 

>co [(7-1 Woo 2 +2J [ 27 /Woo 2 - (7 - 1 )J 


1 

7-i 


( 2 ) 


(3) 


(4) 


(5) 


The ratio p jjp^ depends on vehicle configuration and, for detached shock waves, is not easily or 
accurately computed. The work presented in this paper for the PAET mission is based upon 
experimental evaluations of the pressure coefficient Cp as defined by equation (3). Approximate 
values may be obtained from modified Newtonian theory: Cp = Cp ( cos 2 6, where 

C Pt ~ t J P °°* ~ 1 ^ /^o J p <*)’ and where 0 is the angle between surface normal and vehicle axis. 

Knowledge of local Mach number permits local temperature to be evaluated. The derivation 
follows from equation ( 1 ) and the constancy of total temperature in subsonic flow and across a 
shock „ _ 1 

1 + Moo 2 

T L = — ^ ( 6 ) 

\+ y -^M L 2 

Equations (2) and ( 6 ) can be used to evaluate local velocity and density: 


V T - M .a T =M 


L L 




RTr 


(7) 

( 8 ) 
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Conversion of Kinetic Energy at the Sensor 


Total temperature has two components: the local ambient temperature and the contribution 
of local kinetic energy of the gas stream. Not all the kinetic energy is converted to heat at the 
sensor, and thus even under steady-state conditions a sensor will not record true total temperature. 
This incomplete conversion also implies that the effective temperature for driving a sensor is less 
than total temperature. The effective temperature for driving heat into a sensor— the temperature it 
attains when convective heat transfer goes to zero — is termed recovery temperature. 

The magnitude of recovery temperature is influenced by sensor shape, size, orientation, Mach 
number, flow regime (continuum, slip, or free molecule), and, in some cases, the support-post 
configuration as it affects the flow of gas. Thus, for highest accuracy, a specific sensor should 
ultimately be evaluated in a wind tunnel under conditions representative of those in which it will be 
used. For design purposes, however, data reported in the literature can suffice. A good summary of 
data for cylindrical wires, which rolled butt-welded thermocouples approximate, is given by 
Baldwin et al. (ref. 2), who concluded that for transverse cylinders in continuum and slip flow, 
recovery temperature depends primarily on local Mach number and is almost entirely free of 
Reynolds number dependence for at least 30,000 > Re > 0.10. This range is compatible with the 
Earth and Mars missions and possible sensors. 

For reference, plots of the ratio of 
recovery temperature to total temperature 

RTR = IR (9) 

1 t 

are given in figure 4 as a function of Mach 
number (M<1). These data pertain to 
transverse cylinders in continuum slip flow; 
the data were summarized by Baldwin et al. 

(ref. 2) from refs. 3 and 4. The Mach 
number range 0 to 1 is presented here 
because the flow behind the PAET bow 
wave at the sensor location is subsonic. 



Figure 4. — Values of recovery temperature ratio 
versus Mach number reported in the literature for 
Reynolds numbers greater than 250 and for con- 
tinuum flow. 


Heat Transfer Rate to Sensor 

The rate at which heat is transferred from an air stream to a sensor can be expressed by 

Qc ~ h C (T R ~ r s)C dx (10) 

where h c is the convective heat-transfer coefficient, C dx is surface area of sensor, Tr is local 
recovery temperature, and 7^ is sensor temperature. The convective heat-transfer coefficient is a 
function of sensor shape and size and local flow conditions. Its value is generally presented on a 
dimensionless plot of Nusselt number versus Reynolds number with Mach number as parameter for 
a specific sensor shape. The dimensionless plot of reference 2 for transverse cylinders is reproduced 
in figure 5. 
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Figure 5. — Experimentally determined values of 
Nusselt number versus Reynolds number for trans- 
verse cylinders; curves are from reference 2. 


The convective heat-transfer coeffi- 
cient h c directly influences speed of sensor 
response. The relation between time con- 
stant and h c is 



where r is the time in which a freely 
suspended sensor would respond by 
1 - e ~ l , or 63.2 percent, to a step change 
in recovery temperature; p is sensor 
density; Cp is sensor specific heat; and D is 
sensor diameter. 


Extraneous Heat Sources 

A sensor of the type considered here cannot distinguish among the heating effects of different 
sources. Thus, the effects of solar heating are indistinguishable from those of convective heating 
caused by the gas stream. To ensure that a measurement pertains only to the air stream, effects of 
all other sources of heat should be evaluated to make certain they are negligible. Sources of heat can 
be categorized as: (1) radiant sources, which include the Sun, planet, and vehicle; (2) radiant 
self-emission of the sensor, and (3) conduction through support leads. 

The solar heating effects on a sensor can be expressed 


Q s = oc s F s SCdx ( 12 ) 

where is the absorptance of the sensor for sunlight, F$ is the view factor or effective portion of 
sensor exposed to the flux of one solar constant, S is the solar constant at the planet, and C dx is 
surface area of the sensor. For a sensor exposed broadside to sunlight, F$ equals l/7r; view factors 
are further described in reference 5. Analogously, the heating rates of a sensor due to planetary or 
vehicle infrared radiation are 


Qp = oipFpOTp 4 C dx (13) 

and 

Q v = a v F v oT v *Cdx (14) 

where Tp and T v are the temperatures of planet and vehicle. Sensor self-emission is defined as 
e s oT s 4 C dx, where T s is the temperature and e s is the emittance of the sensor. 

Conduction effects are caused by the relatively massive support posts that do not rapidly track 
total temperature. Thus, heat is conducted from the center of the thermocouple wires to the 
support posts. As a first approximation, conduction effects may be made negligible by means of a 
formula for the steady-state temperature distribution in a wire anchored at both ends (ref. 6). This 
formula states that the temperature error at the center of the wire equals the difference between 
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recovery temperature and support-post temperature divided by cosh m^/2, where L is the length of 
wire and m is defined by 


m =y/h c C/kA (15) 

In this equation, h c is the convective heat-transfer coefficient, C is the circumference of the wire, k 
is thermal conductivity of the wire, and A is the cross-sectional area of the wire. Equation (15) can 
be used in designing the wire sufficiently long to make the temperature error negligible. 


ENERGY EQUATION AND COMPUTER SOLUTION 


The factors that influence measurement were incorporated into an energy equation for relating 
free-stream ambient temperature to sensor response, and the equation was numerically solved by 
computer. To use the resulting computer program for design purposes, one would insert (1) vehicle 
trajectory data, (2) a probable atmospheric temperature profile and thermophysical data, and 
(3) details of sensor design into the program to determine sensor-measured temperature as a 
function of time. Any of the values can be altered to find the effect on sensor response. 


Energy Equation 

The energy equation pertains to a wire * x + AX 

stretched between two support posts, and ; 

is a balance of heat terms on incremental 1 1 

lengths of wire dx as shown in figure 6. The 

equation considers (1) heat transferred Figure 6. — Incremental length of wire on which 

from the air stream, (2) heat conducted energy balance is made 

along the wire, (3) absorbed sunlight, 

(4) absorbed infrared radiation from the planet, (5) absorbed infrared radiation from the entry 
vehicle, (6) sensor-emitted infrared radiation, and (7) changes in heat storage of the sensor. It can 
be written 


(convection + conduction + solar + planet + vehicle)j n p U j- 

- heat storage + (self emission + conduction) ou tp U j- 


That is, 


h c (TR - T S )C dx - kA 


dT { 


+ (a s F s S + 0 i p F p aT p 4 + * v F v oT v 4 )C dx 


ar, 


pCpA -j- dx + eoT s 4 C dx - kA 



(16) 


where T s is sensor temperature, and p, Cp, A, and C are sensor density, specific heat, cross-sectional 
area, and circumference, respectively. Recovery temperature, Tp, is defined by equation (9) and is 
obtained from a plot such as figure 4. The solar, planet, and vehicle radiation terms are defined in 
equations (1 2) through (14). The heat-transfer coefficient h c is obtained from a dimensionless plot, 
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such as figure 5, of Nusselt number versus Reynolds number with Mach number as parameter. 
Equations (2) through (8) provide means of evaluating these dimensionless numbers. Viscosity or 
thermal conductivity are evaluated at total temperature where needed. 

The sensor wire was subdivided into many incremental lengths dx (fig. 6), and the energy 
equation was solved for each length, but subject to defined boundary temperatures at the support 
posts. 


Numerical Solution 


For the numerical solution of the one-dimensional energy equation (16), we first rearrange the 
equation into the form 

d J^ + a(x) d ^+b(x,t,T s )T s = c(t) (17) 

where 


k(x) 

«»■**> 0 


b(x,t,T s ) = 


C[h c {t) + eoT s 3 ] 

pCpA 


C { h c (t)T R (t) + a s F s S + apFpO[ T p (t )] 4 + a v F v o[ T v (t ) ] 4 } 

CW ^pA 


and 


T s T s (x,t ) 

The parabolic partial differential equation (17) is quasilinear, since the coefficient b of the 
nonlinear term is not a function of the derivatives of the dependent variable T s . The solution of this 
equation is uniquely determined in the semi-infinite strip R : (0 <x < L; t > 0} by specifying the 
initial and boundary conditions: 


T s (x,0)=f(x) 0<x<T ) 

T s (0,t) =g 0 (t) t>0 I (18) 

T s (L,t)= gl (t ) t> 0 I 

The analytic solution of the initial boundary-value problem, equations (17) and (18), would 
specify T s and its necessary derivatives as continuous functions of the two independent variables x 
and t. A set of mesh points {(Xj, t n ): xj = j(Ax) where / = 0, 1 , ...,/ + 1 and Ax = L/(J + 1); 
t n =n(At), n = 0, 1, . . .} is introduced on the strip R. The numerical solution is obtained by 
replacing the continuous variables by discrete variables at each mesh point (xj, t n ). 
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The Crank-Nicholson finite-difference scheme is used in finding the numerical solution of 
equations (17) and (18). A description of this scheme, including restrictions for convergence and 
stability, and an error analysis, are given in reference 7. When the Crank-Nicholson procedure is 
applied, functions are evaluated at the point (xp t n+l/2 ), which is halfway between the known 
temperature time level n and the unknown temperature time level (n + 1 ). For subsequent 
notational convenience, let <p(xp t n ) = <f>p n and T s (xp t n ) = Tp n . 


The Crank-Nicholson analogs to the continuous functions are 

T j,n+in 88 2 ( T j,n+i + T j,n) 

(57) “ h " 7 >-») 

(l p) />B+w -2 r ,. B+1 +r / -.,„ +1 ) 

+ (7}-h ,n ~ n + 7/-i ,m)J 


(19) 

( 20 ) 

( 21 ) 


Substituting equations (19) through (21) into equation (17) and rearranging transforms 
equation (17) for the time level (n + 1/2) into a system of algebraic equations 


~ u j^j—i.n+i +v j,n+v 2 ^j,n+i u j^j+i,n+i ~ Sj,n+i /2 1 (22) 


where 


X 


u ) = 2 a i 


X = 


At 

(Ax? 


At 

v j,n+ 1/2 — 1 + Xfly- + “ 2 “ bj n+v2 


and 


Sj,n+ 1/2 u jTj-\ t n + (2 v j,n+v 2 ) + u j^j+i,n + (Af)cy 

The boundary conditions given in equation (18) can be written as 

7o,«+i = £o(?n+i ) 

«>0 (23) 

7 >+i ,72+1 — £i(f/j+i) 

Equations (22) and (23) form a system of equations in which the coefficient matrix is 
tridiagonal. The system is readily solved by using the line-inversion algorithm of reference 8, in 
which the quantities 



di - 


u i 


1 V j,n+V 2 Ujdj—i 


and 


7 = 1 , 2 ,...,/ 


ei = 


C' ^ 

= ^j\n-H/2 + Ujej— X 


1 v j,n+ 1/2 Ujdj— i 


are computed, where d Q = e 0 = 0; 

= ■S'i,n+i/2 + w i ^o, «+i Sj,n+\f2 ~ Sj,n+ 1/2 


/ = 2, 3, . . ., / 


and 


"V.n+i/2 = *$/,«+ 1/2 +u jTj+ \ t n+\ 

The temperature solution for the time level (n + 1) is obtained by backward substitution formulas: 

/ = /- 1 , / - 2 , .... 1 


7>,«+i = e / 

Tj,n+ 1 = djTj+i >n +i + 


, n + 1 


+ 


In figure 7, the centering point for the 
difference equations is indicated by +. The 
circles in the stencil represent the mesh 
points where the temperature values are 
determined. 

The unknown temperature Tj n + V2 
required to calculate n+1/2 in equation 
(22) is initially predicted to be equal to 
Tj n and the solution 7/ ^ +1 is obtained as 
described above. The value 7^+ 1/2 is 
corrected by averaging Tj n and Tj n + t , and 
a new iterated solution 7y jA2+1 is found. This process of iterating over a time step is continued until 
the relative error of two successive solutions is within a prescribed tolerance. 


j-i j j+i 

INCREMENTS ALONG LENGTH OF WIRE 


Figure 7. — Stencil for the Crank-Nicholson finite difference 
scheme. 


The restriction on the time increment (ref. 7) for convergence and stability of the 
Crank-Nicholson finite difference scheme is 

At < ... 

-2a(x) + (Ax) 2 b(x,t,T s ) 

Keller (ref. 7) also shows that for a linear analog the error ej n between the true solution and the 
approximate solution obtained using the Crank-Nicholson method has a second-order bound 

\e jn \<0[(Ax) 2 ] +0[(A/) 2 ] 

The numerical method was coded in Fortran IV(G) using the IBM 360/67 at Ames Research 
Center. 
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PARAMETRIC STUDY 


The computer program was used to evaluate the effects of various parameters on the response 
of a sensor mounted on the PAET or Viking vehicles and flown in proposed Earth and Mars 
trajectories (tables 1 and 2). The method was to vary the value of one parameter at a time and 
observe the effects on the sensor time-temperature history. The effects of (1) sensor diameter, 
(2) sensor length, (3) sensor thermal conductivity, (4) recovery temperature ratio, (5) solar 
radiation, (6) vehicle-emitted infrared radiation, and (7) time of deployment were studied. 

The sensor consisted of a butt-welded chromel-alumel thermocouple wire suspended between 
two support posts (fig. 1). The physical appearance of the weld was that of a continuous piece of 
wire. The wire was 0.0127 cm (0.005 in.) in diameter, and 1.4 cm long. The chromel end of the 
wire was spot-welded to a chromel support post, and the alumel was welded to an alumel post. The 
remainder of the electrical circuit was not specifically considered, but was assumed to be standard 
thermocouple measuring circuit. The sensor was mounted near the outer edge of the forebody, 
S/Rg - 0.985 for PAET and 0.925 for Viking, and was oriented normal to the air flow. It was 
assumed to be housed inside the entry vehicle until a predetermined Mach number was reached, at 
which time it was deployed beyond the boundary layer. 


PAET-Earth Mission 


The Earth-entry computer solutions were initiated, with one exception, when the entry vehicle 
had slowed to a Mach number of 2. The sensor on PAET was to be deployed at this Mach number. 
The values assigned to the various parameters in the computer solutions are given in table 3. In 
addition to these, (1) recovery temperature ratios were obtained from Scadron and Warshawsky’s 
data (ref. 4), plotted in figure 4; (2) heat-transfer coefficients were obtained from figure 5; 
(3) pressure coefficients, Cp were obtained from experimental data, figure 8; (4) support-post 
temperatures were obtained from a computer solution of freely supported posts at flow conditions 
corresponding to the sensor location, figure 9; and (5) surface temperatures of the vehicle forebody, 
which radiates to the sensor, were obtained from predictions, figure 10. 



Figure 8. — Pressure coefficient versus Mach number for Figure 9. — Support-post temperatures for the 

the PAET vehicle at zero angle of attack; sensor PAET-Earth study, 

located at S/Rg = 0.985. 
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Figure 10. — Forebody temperature, near the sensor, 
versus time. 


TIME, sec 



(a) Effect of sensor diameter. 

TIME, sec 



(b) Effect of sensor length and thermal conductivity. 


Figure 11. — Response of sensors deployed at = 2 
during PAET-Earth mission. 


The time response of the 
0.0127-cm-dia sensor to be flown on PAET 
is illustrated in figure 1 1 (a). The response 
curve corresponds to computer run 1 in 
table 3, and represents the nominally ex- 
pected values of the various parameters. 
T his figure is a plot of sensor temperature 
versus time, with time beginning at the 
instant of deployment and continuing to 
30 sec. This time span is considered illustra- 
tive of the problems of measurement. In 
actual flight, measurements will be made 
over a longer time. The dashed curve 
represents recovery temperature, which is 
the temperature directly available for 
driving convective heat into the sensor. 

The sensor temperature is seen to 
increase rapidly after deployment, where 
its initial temperature is 300° K, and in 0.4 
sec reaches recovery temperature. There- 
after, however, the sensor lags this tempera- 
ture because of its own thermal mass and 
because vehicle deceleration causes the 
recovery temperature to fall rapidly. At 
t = 1 sec, the difference between sensor and 
recovery temperatures is 6° K; at t = 6 sec, 
3° K; and at t = 30 sec, less than 1° K. Thus 
the maximum difference occurs at about 1 
sec, and thereafter continually decreases 
with time. The useful measurement period 
begins at t = 0.4 sec, when Moo - 1 -92 and 
the altitude is 30.2 km. 

Figure 11(a) also illustrates the effect 
of sensor diameter on response. The curve 
labeled D/5 is for a 0.00254-cm 
(O.OOl-in.)-dia sensor, and the one labeled 
2D is for a 0.0254-cm (O.OlO-in.)-dia 
sensor. These correspond to computer runs 
2 and 3 in table 3. The smaller sensor has a 
decided advantage in speed of response and 
overall accuracy. It reaches recovery tem- 
perature in about 0. 1 sec, and thereafter its 
temperature difference with respect to 
recovery temperature is less than 1.5° K. 
The larger sensor has a very slow time 
response, requiring almost 0.8 sec to reach 
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recovery temperature. It has a maximum 
difference of 10° K, and this only slowly 
decreases with time. The reasons the inter- 
mediate size sensor was selected for the 
PAET mission are that its time response 
was adequate for objectives, and that it was 
considerably more rugged, for handling 
purposes, than the smaller one. 

The effects of varying the length and 
thermal conductivity of the sensor are 
illustrated in figure 1 1(b). Three lengths of 
wire — corresponding to L, L/2, and 2 L — 
and the case of zero thermal conductivity 
were considered. The purpose of these 
studies was to evaluate how closely the var- 
ious lengths of wire correspond to the ideal 
situation of an infinitely long wire, or zero 
thermal conductivity. Thus, in this figure, 
the zero thermal conductivity case is con- 
sidered ideal. At about 0.4 sec after deploy- 
ment, the difference between the ideal case 
and the other sensors is a maximum. The 
maximum difference is about 8° K and 
depends on sensor length. The reason for 
these large temperature differences is that 
the wire rapidly responds to the high re- 
covery temperature, but the more massive 
support posts do not respond so fast, with 
the consequent loss of heat by conduction 
to the support posts. Fortuitously, how- 
ever, the L/2 sensor would fall closer to 
true recovery temperature for relatively 
early times. At later times, the temperature 
of the L - 1.4-cm sensor is effectively that 
of an infinitely long sensor wire. 

The effects of solar, vehicle-emitted, 
and sensor-emitted radiation on the 
accuracy of measurement are illustrated in 
figure 11(c). Each radiation level was 
assigned a value of zero, corresponding to 
computer runs 7, 8, and 9 in table 3, and 
these were compared with the response 
curve for the nominally expected radiation 
levels (computer run 1 ). The responses for 
zero values of solar and sensor-emitted 
radiation were so nearly identical to those 


TIME, sec 



(c) Effect of solar, vehicle-emitted, and 
sensor-emitted radiation. 


TIME, sec 


6 10 14 18 22 26 30 , 



(d) Effect of recovery temperature ratio. 


Figure 11. — Concluded. 



Figure 12. — Local Mach number, at sensor location, 
versus time; for PAET-Earth mission. 
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(a) Effect of sensor diameter. 



(b) Effect of sensor length and thermal conductivity. 

Figure 13. — Response of sensors deployed at Moo = 3 
during PAET-Earth mission. 


for the nominally expected levels that for 
the most part these curves are shown as 
one. The relation of one to the other is 
shown slightly exploded at about 0.4 sec 
where the small differences start to appear 
(the actual total spread of the four curves 
at t = 0.4 sec was found to be only 1 .6° K). 
The largest error was for the case of zero 
vehicle radiation, but even this error was 
less than about 1° K. It can be concluded 
that for the PAET mission, and deploy- 
ment at Mach number 2, the influence of 
radiation on measurement is negligible. 

Figure 11(d) illustrates the effects of 
varying recovery temperature ratio. The 
cases of RTR = 1 (computer run 10) were 
compared with those for RTR from figure 
4, ref. 4 (computer run 1). The effect is 
quite large, amounting to the order of 
10° K, and only slowly decreases with 
time. The reason for this behavior can be 
deduced from figures 4 and 12. Figure 12 
shows that local Mach number (at the 
sensor location behind the shock) is fairly 
constant but does slowly decrease with 
time. When the free-stream Mach number is 
1 or 2, the local Mach number is only 
about 0.65. It can be seen from figure 4, 
that at Mach 0.65, RTR ~ 0.975. Thus, the 
difference between the two response curves 
at early times after deployment amounts to 
about 0.025 of total temperature, or 
approximately 10° K. Since local Mach 
number decreases only slowly with time, 
the temperature error too will decrease 
only slowly. Consequently, it is important 
to accurately know the value of the 
recovery temperature ratio. 

The corresponding computer runs for 
deployment of the sensor at Mach number 
3 during the PAET-Earth mission are 
illustrated in figures 13(a)-(d), which have a 
time scale going to 4 sec after deployment. 
Note that the Mach number decreases to 2 
at about t = 3 sec. These figures illustrate 
that error of measurement in the first 3 sec 
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would be twice that for the Mach 2 
deployment. The basic reason is that total 
and recovery temperatures are higher, and 
changing much more rapidly than at Mach 
number 2 (fig. 3), which influences the 
ability of the sensor to track the change. It 
can be seen that the difference between 
recovery temperature and sensor tempera- 
ture for the 0.0127-cm wire, at a time of 
1 sec, is about 12° K. For the 2D case, the 
error approaches 22° K. The influence of 
sensor length is more critical than before 
because of the increased temperature 
difference between the center of the sensor 
wire and support posts. Errors of the order 
7° K result. The effects of radiation are 
again small, but larger than for Mach 2. The 
corresponding comparison of recovery 
temperature ratios produces an error of 
about 0.025 total temperature, or about 
14° K. Thus the tracking ability of the 
sensor is degraded above Mach number 2. 



(c) Effect of solar, vehicle-emitted, and 
sensor-emitted radiation. 


Viking-Mars Mission 

An abbreviated form of the computer 
program was used for the Viking-Mars 
study. First, local Reynolds numbers at the 
sensor location were input to the computer 
(fig. 14); these numbers were used because 
experimental values of pressure coefficient 
Cp were not available, but Reynolds 
numbers were available from a flow-field 
study. Second, Mach number effects on 
heat-transfer coefficients (fig. 5) were 
neglected, and all values were computed at 
an arbitrary Mach number of 0.65. (It will 
be recalled that flow at the sensor location 
is always subsonic.) Third, Mach number 
effects on recovery temperature ratios (fig. 
4) were neglected, and all ratios were 
arbitrarily taken as 1.0. Fourth, the 
temperatures of the support posts were 
assumed to be constant. The values 
assigned to the various parameters are listed 
in table 4; the temperature of the vehicle 



(d) Effect of recovery temperature ratio. 


Figure 13. — Concluded. 
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Figure 14. - Reynolds number versus time, at sensor 
location for Viking-Mars mission. 



(a) Effect of sensor diameter. 



(b) Effect of sensor length and thermal conductivity. 

Figure 1 5. - Response of sensors deployed at Moo = 3 
during Viking-Mars mission. 


forebody, which radiates to the sensor, is 
illustrated in figure 10. The sensor 
Reynolds numbers, local Mach number, 
heat-shield temperatures, and other param- 
eters of the flow for these calculations 
were provided by Mr. H. C. Norman of 
Martin Marietta Corporation from results 
obtained with a flow-field computer 
program. 

The results of the parametric study 
are shown in figures 15(a)-(c) for deploy- 
ment at Mach number 3. A 6-sec time span 
was selected as illustrative of the measure- 
ment problems. The response of a PAET 
sensor deployed during the Viking-Mars 
mission, is shown in figure 15(a); 
D = 0.0127 cm. It reaches the recovery 
temperature in about 1.5 sec and thereafter 
tracks it quite well, differing by only about 
1° K. Conduction effects for the selected 
support-post temperature of 300° K are 
just sufficient to counterbalance the lag in 
the sensor temperature. The response of 
sensors having D/5 and 2D diameters is also 
shown. The D/5 sensor responds in less 
than 0.1 sec and thereafter tracks recovery 
temperature within 1° K. The 2D sensor, 
however, responds very slowly and after 
3 sec has errors of the order of 1 2° K. 

Figure 15(b) shows that the sensor 
selected for the PAET mission is not quite 
long enough to eliminate conduction 
effects for the Viking-Mars mission; if the 
sensor were shortened by a factor of 2, 
temperature errors of the order 45° K 
would result. Because of conduction 
effects, the sensor temperature will be 
influenced by the specific values of 
support-post temperatures. 

The effect of radiation on measure- 
ment accuracy is shown in figure 15(c). 
Solar radiation has nominally no effect due 
to the long distance of Mars from the Sun. 
Sensor self-emission and vehicle forebody 
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(c) Effect of solar, vehicle-emitted, and sensor- 
emitted radiation. 

Figure 15. — Concluded. 


radiation each produce errors of the order 1° K, and in opposite directions, so these effects tend to 
cancel one another. 


Deployment at Mach number 4 is illustrated in figures 16(a)~(c). The PAET sensor and the D/5 
diameter sensor track recovery temperature quite well, but the 2D sensor has errors of the order 
16° K. The effect of sensor length is very prominent, and the L/2 sensor would produce errors of the 
order 60° K. Figure 1 6(c) shows that the effect of sensor self-emission is about 4° K. 




(b) Effect of sensor length and thermal conductivity. 


Figure 16. — Response of sensors deployed at Moo = 4 during Viking-Mars mission. 
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(c) Effect of solar, vehicle-emitted, and sensor-emitted radiation. 
Figure 16. — Concluded. 


CONCLUSIONS 


An analytical study, such as the one described here, must be conducted to define the behavior 
and evaluate the accuracy of an atmospheric recovery-temperature sensor to be used for specific 
atmospheric-measurement missions. The following mission-specific factors, shown to influence the 
sensor design and response were included in the analysis: the relation of free-stream conditions to 
local Mach and Reynolds numbers at the sensor location, the heat-transfer characteristics between 
the sensor and the atmosphere, and extraneous heat sources such as conduction and radiation. The 
computer program developed for this study should be valuable for analyzing the behavior of sensors 
for future missions. 

For the sensor selected to be flown on PAET and deployed at Mach number 2, the potential 
source of greatest error is an inaccurate knowledge of the recovery temperature ratio. The effects of 
radiation are negligible, and the sensor was designed to minimize conduction effects. For 
deployment at a higher Mach number (3), the sensor does not track recovery temperature as well 
because the temperature is higher and changes more rapidly with time. At these higher 
temperatures, errors in the recovery temperature ratio also produce larger errors in the deduced 
ambient temperature. Deployment at Mach number 3 approximately doubles the measurement 
error (during the increased time of measurement). 

If a sensor similar to that designed for PAET, but having only half the wire length, were used 
for the Viking-Mars mission, the resulting errors due just to conduction effects would be of the 
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order 45° to 60° K. Thus, one cannot assume that a sensor proved adequate for one mission will be 
adequate for another. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, California 94035, April 10, 1972 
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TABLE 1 - TRAJECTORY AND PLANETARY DATA FOR PAET-EARTH MISSION 


Time, 

sec 

Altitude, 

km 

Ambient 

temperature, 

°K 

. Ambient 
pressure, 
N/m 1 2 

Ambient 

density, 

kg/m 3 

Mach 

number 

Total 

temperature, 

°K 

Ratio of 
specific heats 

0 

31.6 

228 

942 

0.0144 

3.00 

639 

1.4 

2.5 

30.6 

227 

1100 

.0169 

2.17 

441 


5 

29.8 

226 

1241 

.0191 

1.66 

351 


10 

28.5 

225 

1497 

.0232 

1.10 

280 


20 

26.5 

223 

2017 

.0315 

.769 

249 


30 

24.7 

221 

2668 

.0420 

.649 

240 


40 

23.0 

220 

3475 

.0551 

.569 

234 


50 

21.4 

218 

4429 

.0708 

.505 

229 


100 

15.7 

217 

10914 

.175 

.309 

221 


150 

11.8 

217 

20035 

.322 

.225 

219 


200 

8.82 

231 

31634 

All 

.178 

232 


250 

6.31 

247 

45267 

.638 

.149 

248 


300 

4.11 

261 

60806 

.810 

.128 

262 


350 

2.14 

274 

78151 

.993 

.113 

275 


400 

.347 

286 

97225 

1.18 

.101 

286 


410 

.007 

288 

101242 

1.22 

.099 

288 

- - . __! 



TABLE 2 - TRAJECTORY 1 AND PLANETARY DATA FOR VIKING-MARS MISSION 


Time, 

sec 

Altitude, 

km 

Ambient 

temperature, 

°K 

Ambient 

pressure, 

N/m 2 

Ambient 

density, 

kg/m 3 

Mach 

number 

Total 

temperature, 

°K 

Ratio of 
specific heats 

0 

23.8 

197 

113 

0.00280 

3.98 

790 

1.38 

5 

22.8 

198 

125 

.00308 

3.47 

651 

1.38 

10 

21.7 

199 

136 

.00334 

3.05 

550 

1.38 

20 

19.8 

200 

163 

.00397 

2.39 

419 

1.38 

30 

17.9 

202 

193 

.00465 

1.93 

345 

1.38 

40 

16.1 

204 

225 

.00537 

1.60 

302 

1.38 

50 

14.4 

206 

265 

.00626 

1.35 

275 

1.37 

100 

5.4 

213 

572 

.0131 

.85 

242 

1.37 

130 

0 

230 

905 

.0190 

.72 

252 

1.36 


1 Ballistic entry with entry velocity of 4,570 m/sec (15,200 ft/ sec) ballistic coefficient of 47 kg/m 2 

(0.3 slug/ ft 2 ), and an entry path angle 16° below the horizontal. 
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TABLE 3 - VALUES OF PARAMETERS USED IN COMPUTER PROGRAM FOR PAET-EARTH STUDY 


Computer run number 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Sensor 











Diameter, cm 

0.0127 

Dl 5 

2D 

0.0127 







Length, cm 

1.4 



LI 2 

2L 

1.4 





Emittance 

0.14 








0 

0.14 

Absorptance for sunlight 

0.40 






0 

0.40 



Absorptance for planet radiation 

0.14 










Absorptance for vehicle radiation 

0.20 







0 

0.20 


View factor for sunlight 

1/tt 










View factor for planet radiation 

0.50 










View factor for vehicle radiation 

0.50 










Specific heat, J/kg°K 

460 










Density, kg/m 3 

8666 





0 





Thermal conductivity, chromel, W/m°K 

19.0 





0 

19.0 




Thermal conductivity, alumel, W/m°K 

29.4 






29.4 




Support post temperature 

See fig. 9 










Recovery temperature ratio 
Planet 

See fig. 4 









1.0 

Temperature, °K 

273 










Distance from Sun, AU 

1.0 













TABLE 4 - VALUES OF PARAMETERS USED IN COMPUTER PROGRAM FOR VIKING-MARS 

STUDY 


Computer run number 

1 

2 

3 

4 

5 

6 

7 

8 

9 

Sensor 










Diameter, cm 

0.0127 

D/5 

2D 

0.0127 






Length, cm 

1.4 



L/ 2 

2 L 

1.4 




Emittance 

0.14 








0 

Absorptance for sunlight 

0.40 






0 

0.40 


Absorptance for planet radiation 

0.14 









Absorptance for vehicle radiation 

0.20 







0 

0.20 

View factor for sunlight 

1/tt 









View factor for planet radiation 

0.50 









View factor for vehicle radiation 

0.50 









Specific heat, J/kg °K 

460 









Density, kg/m 3 

8666 









Thermal conductivity, chromel, W/m° 

19.0 





0 

19.0 



Thermal conductivity, alumel, W/m U K 

29.4 





0 

29.4 



Support post temperature, °K 

300 









Recovery temperature ratio 

1.0 









Planet 










Temperature, °K 

280 









Distance from Sun, AU 

1.52 
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